More on the infrared renormalization group limit cycle in QCD 

E. EpelbauiiE 

Jefferson Laboratory, Theory Division, Newport News, VA 23606, USA 

H.-W. HammeJll 

Helmholtz-Institut fiir Strahlen- und Kernphysik (Theorie), 
Universitdt Bonn, Nufiallee 14-16, D-53115 Bonn, Germany 

Ulf-G. MeiBneJl 

Helmholtz-Institut fiir Strahlen- und Kernphysik (Theorie), 
Universitdt Bonn, Nufiallee 14-16, D-53115 Bonn, Germany and 

Institut fiir Kernphysik (Theorie), Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

A. Noggci 

Institut fiir Kernphysik (Theorie), Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

(Dated: February 2, 2008) 

Abstract 

We present a detailed study of the recently conjectured infrared renormalization group limit 
cycle in QCD using chiral effective field theory. It was conjectured that small increases in the 
up and down quark masses can move QCD to the critical trajectory for an infrared limit cycle in 
the three-nucleon system. At the critical quark masses, the binding energies of the deuteron and 
its spin-singlet partner are tuned to zero and the triton has infinitely many excited states with 
an accumulation point at the three-nucleon threshold. We exemplify three parameter sets where 
this effect occurs at next-to-leading order in the chiral counting. For one of them, we study the 
structure of the three-nucleon system in detail using both chiral and contact effective field theories. 
Furthermore, we investigate the matching of the chiral and contact theories in the critical region 
and calculate the influence of the limit cycle on three-nucleon scattering observables. 
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I. INTRODUCTION 



The renormalization group (RG) is an important tool in many areas of physics. Its apph- 
cations range from critical phenomena in condensed matter physics to the nonperturbative 
formulation of quantum field theories for elementary particles P. The requirement of the 
invariance of low-energy observables under changes of the ultraviolet cutoff A generates a RG 
flow on the multidimensional space of coupling constants g for operators in the Lagrangian: 

Aj|g = /3(g). (1) 

Depending on the physical application, the RG flow can be considered in the infrared (A 0) 
or in the ultraviolet (A —>■ oo) limits. The simplest solution to the RG equations is a 
fixed point g* which satisfies /3{g*) = 0. An important signature of an RG fixed point is 
scale invariance: symmetry with respect to the coordinate transformation r — > Ar for any 
positive number A. This symmetry implies that dimensionless variables scale as powers of the 
momentum scale. Scale-invariant behavior at long distances, as in critical phenomena, can be 
explained by RG fiow to an infrared fixed point. Scale- invariant behavior at short distances, 
as in asymptotically-free field theories, can be explained by RG fiow to an ultraviolet fixed 
point. 

However, RG equations can also exhibit more complicated solutions. The possibility of 
RG fiow to a limit cycle was proposed by Wilson in 1971 A limit cycle is a 1-parameter 
family of coupling constants g*(^) that is closed under the RG fiow and can be parametrized 
by an angle < ^ < 27r. The RG fiow carries the system around a complete orbit of the limit 
cycle every time the ultraviolet cutoff A increases by some factor Aq. One of the signatures 
of an RG limit cycle is discrete scale invariance: symmetry with respect to the coordinate 
transformation r — > Agr only for integer values of n. 

Before the fundamental theory of the strong interactions was known, Wilson had sug- 
gested that limit cycles might be relevant to the high-energy behavior of the strong in- 
teractions of elementary particles [2]. However, Quantum Chromodynamics (QCD) and 
asymptotic freedom were discovered soon thereafter Q and the high-energy behavior of the 
strong interactions was explained by an ultraviolet fixed point in QCD. The low-energy 
structure of the strong interactions, however, is much more complicated and not dominated 
by the fixed point. 

The low-energy sector can be described by exploiting the approximate chiral symmetry 
of QCD using chiral effective field theory (EFT) methods Chiral EFT is a powerful 

tool for analyzing the properties of hadronic systems at low energies in a systematic and 
model-independent way. It is formulated in an expansion around the chiral limit of QCD 
which governs low-energy hadron structure and dynamics. In chiral EFT's the quark mass 
dependence of the operators in the effective Lagrangian is included explicitly. Over the past 
15 years, considerable progress has been made in understanding the structure of the nuclear 
force in this framework 



The quark mass dependence of the chiral nucleon-nucleon (NN) interaction was studied 
with the primary aim to understand the chiral limit of nuclear physics 0, ^ 3| . While this 



question has not been fully resolved, these studies found that the inverse scattering lengths 
in the ^Si-^Di and ^Sq channels both vanish if one extrapolates away from the physical limit 
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to slightly larger quark masses.^ Subsequently, it was pointed out that QCD is close to the 
critical trajectory for an infrared RG limit cycle in the 3-nucleon sector. This led to the 
conjecture, that QCD could be tuned to the critical tra ject ory by small changes in the up 
and down quark masses away from their physical values [15|. The proximity of the physical 
quark masses to these critical values explains the successful description of the low-energ 
three-nucleon problem in terms of zero-range forces between nucleons initiated long ago [if 
The effective-field-theory formulation of this program exhibits an ultraviolet RG limit cycle 
[l7| . The proximit y o f physical QCD to the critical trajectory implies that the ultraviolet 
limit cycle of Ref. [l7| is not just an artifact of the EFT but hints towards a limit cycle in 
QCD. 

The connection between the limit cycle and three-nucleon observables is established by 
the Efimov effect [18] which occurs in the three-body sector of nonrelativistic particles with 
a resonant short-range S-wave two-body interaction. The strength of the interaction is 
governed by the S-wave scattering length a. If a is large and positive, there is a shallow 
two-body bound state with binding momentum k = 1/a, if a is large and negative there is 
a shallow virtual state characterized by the momentum scale l/|a|. Efimov showed that if 
|a| is much larger than the range tq of the interaction, there are shallow three-body bound 
states whose number increases logarithmically with |a|/ro. In the resonant limit a — ±cxd, 
there are infinitely many shallow three-body bound states with an accumulation point at the 
three-body scattering threshold. If the particles are identical bosons, the ratio of the binding 
energies of successive states rapidly approaches the universal constant Ag ~ 515. Efimov 
also showed that low-energy three-body observables for different values of a are related by 
a discrete scaling transformation in which a AqO, where n is an integer, and lengths and 
energies are scaled by the appropriate powers of Aq 3 3- The mathematical connection 
between the Efimov effect and RG limit cycles was first pointed out in Ref. [20]. 

The Efimov effect can also occur for fermions with at least three distinct spin or isospin 
states and therefore applies to nucleons as well. The spin-singlet and spin-triplet np scat- 
tering lengths are oisq = —23.8 fm and aas^ = 5.4 fm. They are both significantly larger 
than the effective range, which is tq = 1.8 fm in the spin-triplet channel. Efimov used 
this observation as the basis for a qualitative approach to the three-nucleon problem [21j . 
A convenient implementation of this program is given by the so-called pionless or contact 
EFT j^, [i^l • Nucleons are described as point particles with zero-range interactions whose 
strengths are adjusted to reproduce the scattering lengths asg^ and aig^. The effective range 
and higher order terms in the low-energy expansions of the phase shifts are treated as per- 
turbations. This approach works well in nuclear few-body systems dominated by momenta 
small compared to M^-. In the triton channel, the Efimov effect makes it necessary to include 
a three-body force at leading order in the power counting [l^. The three-body force can 
be fixed by using one three-body datum as input. All other three-body observables can 
then be predicted. The structure of this EFT is much simpler than the chiral EFT and the 
computational effort is considerably smaller. Since the contact EFT is based on an expan- 
sion around the limit of infinite scattering length, it is particularly well suited to describe 
processes governed by the large scattering length. As a consequence, the chiral and contact 
EFT's can mutually complement each other. A first exploratory study of the infrared limit 



Due to the nuclear tensor force, the '^Si and '^Di channels are coupled. This mixing is included in the 
chiral EFT calculations while it appears as a higher order effect in the contact EFT discussed below. For 
simplicity, we will only refer to the '^Si and ^Sq partial waves in the following. 
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cycle in QCD in the contact EFT was carried out in Ref. 3|- The quark mass dependence 



of the nucleon-nucleon scattering lengths from Ref. fli*] was used as input in this calculation. 

In this paper, we study the possibility of an infrared limit cycle in QCD in more depth. 
We use both chiral EFT and pionless EFT and combine the strengths of both approaches. 
We use the chiral EFT to calculate the bound state spectrum of the triton in the vicinity of 
the limit cycle and study how well it is approximated by the contact EFT. Furthermore, we 
calculate various three-body scattering observables in pionless EFT and illustrate how they 
are affected by the limit cycle. 

II. CHIRAL EFFECTIVE FIELD THEORY 

The quark mass dependence of the chiral A^A^ interaction was calculated to next-to- 
leading order (NLO) in the chiral counting in Refs. jl^. [isL At this order, the quark 



mass dependence is synonymous to the pion mass dependence because of the Gell-Mann- 
Oakes-Renner relation: 

Ml = -{mu + md){^\uu\^)/Fl, (2) 

where (OImmIO) ~ (—225 MeV)^ is the quark condensate. In the following, we will therefore 
refer only to the pion mass dependence which is more convenient for nuclear applications 
and treat the pion mass as a parameter that can be varied by adjusting the values of the 
quark masses. In the work of Refs. 0, 0, it was found that the scattering lengths in 



the ^Si-^Di and ^Sq channels both vanish in the pion mass region around 200 MeV. 

In this work, we study the structure of the nuclear three-body system in this region of 
pion masses around 200 MeV. We will use the chiral A^A^ potential constructed from EFT 
using the method of unitary transformations [14]. To next-to-leading order (NLO) in the 
chiral power counting, this potential can be written as: 

y^Lo = ^^'''^ + V^'''' + V^""^' , (3) 

where V^^^, V"^^^, and V'^°'^^ refer to the one-pion exchange, two-pion exchange, and contact 
potentials, respectively. They are given by the expressions: 

^OP. . 4|| + - ^..1 ■ ^^^^m^ , (4) 



AMl{^g\ - Ag\ - 1) + q^(23^1 - 10^1 - 1) + 



4M2 + q2_ 



+ M^i-q^2-q-q f^i-f72^ , (5) 
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= Cs + Cria, ■ a,) + ( Ds - .^^A^FICt - hg\ + 2) In ■ 



+ ^C5^^^ ■ (k X q) + C6(q • ai)(q ■ ^2) + C7(k ■ ai)(k ■ ^2) , (6) 

with gA and the physical values of the nucleon axial coupling and pion decay constant, 
resp ectively and Jig a low-energy constant related to the Goldberger-Treiman discrepancy 



14| . The symbols (jj (tj), z = 1,2 indicate the spin (isospin) operators for particle i 
and the Ci^ . j, Cs;r-, and Ds,t are low-energy constants (LEC's) to be determined from 
fits to nucleon-nucleon data. Further, q denotes the momentum transfer of the nucleon, 
i.e. q = p' — p, where p' and p are final and initial nucleon momenta, while k = (p' -(- p)/2. 
Here and in what follows we denote the value of the variable pion mass by M^^ in order to 
distinguish it from the physical value denoted by M'^^^ = 139.6 MeV. Furthermore, 

^'^'^ |q| 2M^ ' ^ ^ 

and A represents the relative shift in the ratio (^aZ-^tt compared to its physical value: 

{9a/ Fn)M^ - {9a/ F^) MS'^y 



A = 



{9a/F^)j., 



phys 



where the low-energy constants diQ and Z4 are defined as in Ref . |1^ . Note that in the TPEP, 
we only take into account the explicit M^-dependence and use the physical values for qa 
and Ftt. This is sufficient at NLO since any shift in qa and F-^ for a different value of M^r 
in the TPE is a N^LO effect. We also incorporate the leading isospin-breaking corrections 
due to the pion mass difference in the one-pion exchange potential [2^ but do not consider 
an independent variation of rriu — rrid. 

The constants Cst and Dst are related to the Cst used in 0] via 



Cs,T = Cs,T + {Mfynms,T. (9) 

Note further that the short-range terms of the type In M^r in Eq. (jHl) result from the 
two-pion exchange as well as from the renormalization of the leading-order contact forces 
by pion loops. It is important to stress that renormalization of the LECs Cs, Ct, Ci,...7 
due to pion loops does not depend on the pion mass and thus is of no relevance for this 
work. The potential V{p', p) is multiplied by the regulating functions /r(|p |), /r(|p'|) in 
order to cut off the large momentum components in the Lippmann-Schwinger equation. In 
this study, we use the same exponential function /ij(|p |) = exp[— p'^/A^] as in Ref. ji^ and 
restrict ourselves to the cut-off A = 540 MeV^. 



^ The cutoff dependence is mild in the very low energy regime, in which we are interested here. For a 
discussion of larger cutoffs, see Ref. 0- 
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In principle, these equations determine the pion mass dependence of the chiral A^A^ 
potential uniquely. However, the extrapolation away from the physical pion mass generates 
errors. The dominating source are the constants Cs,t and Ds^t which cannot be determined 
independently from fits to data at the physical pion mass. A smaller effect is due to the 
error in the LEG di^ which is enhanced in Eq. (jSI) as one moves away from the physical 
pion mass. Both effects generate increasing uncertainties as one extrapolates away from the 
physical point. Note also that we do not include explicit A (1232) degrees of freedom. It 
would be interesting to see if and how our results would be modified in a theory with explicit 
A's 0. 

Following Ref. l3|, we use di% = —1.23 GeV~^ which is the average of three values given 



in Ref. j27| (See also Ref. [28^). The LEG dig = -0.97 GeV'^ is determined from the ob- 
served value of the Goldberger-Treiman discrepancy while the two remaining constants Ds 
and Dt are unknown. The size of these two constants can be constrained from natural- 
ness arguments. In Ref. 3|, it was argued that the corresponding dimensionless constants 



F^A'^Ds^T can be expected to satisfy the bounds: 

-3 < FlAlDs,T < 3 , (10) 
where ~ 1 GeV is the chiral symmetry breaking scale. We note that Refs. jl^ . 13l | do allow 



for a larger variation of these LEG's. The LEG's that we are going to choose below are within 
our more restrictive range. Therefore, the discussion on the appropriate scaling of the Ds,t 
is not relevant for this study. Furthermore, these bounds are in agreement with resonance 
saturation estimates and similar conditions are obeyed by the known constants |29|| . For the 
constants Cs^Ti e.g., we find Cs = —120.8 GeV~^ and Ct = 1.8 GeV~^ corresponding to the 
dimensionless coefficients F'^Cs = —1.03 and F^Ct = 0.02, respectively. The unnaturally 
small value of F'^Ct is a consequence of the approximate Wigner SU(4) symmetry. (For a 
discussion of this issue in the pionless EFT, see Ref. j|40|]). 

The ranges from Eq. ()10|) were used to estimate the extrapolation errors of two-nucleon 
observables like the deuteron binding energy and the spin-singlet and spin-triplet scattering 



lengths in Ref. [IJ]. In the exploratory study of the three-nucleon system [15[, the mean 
values of these error bands were used as input for the three-body calculations in the contact 
EFT. Even though both scattering lengths were large for the mean values, they did not 
become infinite at the same value of the pion mass and there was no exact limit cycle for 
this choice of parameters. 

Here we take a different approach and search for sets of values for Ds and Dt that lie 
within the bound given by Eq. (fTUj) and cause the spin-singlet and spin-triplet scattering 
lengths to become infinite at the same value of the pion mass. For this purpose, it is more 
convenient to use the partial wave projected constants 

Das, = 47r {Ds + Dt) , 

Diso = 47r [Ds - 3Dt) . (11) 
The dimensionless constants and asg, defined as 

ais, = F^AlDisJilQir) and mg, = F^AlD.sJiSir) , (12) 
then satisfy the same bound as in Eq. (fTI]|) . i.e.: 

— 3 < «3s, < 3 , 

-3 < aiso < 3 . (13) 
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FIG. 1: Inverse of the S-wave scattering lengths in the spin-triplet and spin-singlet nucleon-nucleon 
channels for solution (a) as a function of the pion mass M^. The vertical dotted line indicates the 
critical pion mass M^****. 

The values of the parameters ass^ and aisg can be chosen independently anywhere in the 
above intervals without violating naturalness and without affecting physics at the physical 
value of the pion mass. 

The values of the pion mass where the inverse scattering lengths in the spin-triplet and 
spin-singlet channels vanish simultaneously are called critical pion masses M^^**. At NLO in 
the chiral counting, it is possible to find parameter sets with critical pion masses in the range 
175 MeV ^ M^''** < 205 MeV. For example, for the three exemplifying values aa^^ = ±2.5 
and 0.0, we obtain the following critical parameter sets: 

(a) ass, = -2.5 and aig, = 2.138598 =^ M^"* = 197.8577 MeV, 

(b) ass, = 0.0 and aig,, = 1.955709 =^ M^"* = 186.3276 MeV, 

(c) ass, = 2.5 and aig, = 1.776665 =^ M^"* = 179.0417 MeV. 

We note that it is unlikely that physical QCD will correspond to any of the solutions (a)-(c). 
However, in Ref. 3| it was conjectured that one should be able to reach the critical point by 



varying the up- and down-quark masses m„ and independently because the spin-triplet 
and spin-singlet channels have different isospin. A more detailed investigation is needed in 
order to test this conjecture. However, many aspects of the limit cycle are universal and 
do not depend on the exact parameter values [i^. Therefore, we study the structure of the 
three-nucleon system near the critical pion mass for solution (a) in the remainder of this 
paper in more detail. The universal aspects of the three-body observables do not depend on 
the details of the solution we choose. 

The inverse scattering lengths in the spin-triplet and spin-singlet channels in the vicinity 
of the limit cycle for solution (a) are shown in Fig. ^ As promised, the inverse scattering 
lengths vanish at the critical value of the pion mass M^"* = 197.8577 MeV. For pion masses 
below the critical value, the spin-triplet scattering length is positive and the deuteron is 
bound. As the inverse spin-triplet scattering length decreases, the deuteron becomes more 
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FIG. 2: Binding energies of the triton ground state and the first two excited states as function of 
the pion mass M,r- The circles (ground state), squares (first excited state), and diamonds (second 
excited state) give the chiral EFT result, while the dashed lines are calculations in the pionless 
theory. The vertical dotted line indicates the critical pion mass M^"^. The thresholds for the 
three-body states are given by the solid lines. 

and more shallow and finally becomes unbound at the critical mass. Above the critical pion 
mass the deuteron exists as a shallow virtual state. In the spin-singlet channel, the situation 
is reversed: the "spin-singlet deuteron" is a virtual state below the critical pion mass and 
becomes bound above. The pion mass dependence of the two scattering lengths shown in 
Fig. ^ will be used as input for the calculations in the contact EFT in the next section. 

From the solution of the Faddeev equations with solution (a) for the A^A^ potential, we 
obtain the binding energies of the triton and the first two excited states in the vicinity of 
the limit cycle (See Ref. [s^ for details). The binding energies are given in Fig. |21by the 
circles (ground state), squares (first excited state), and diamonds (second excited state). The 
solid lines indicate the neutron-deuteron {M^, < M^"*) and neutron-spin-singlet-deuteron 
(Mtt > M^^**) thresholds where the three-body states become unstable. Directly at the 
critical mass, these thresholds coincide with the three-body threshold and the triton has 
infinitely many excited states. The dashed lines are calculations in the pionless theory and 
will be discussed in detail below. The binding energy of the triton ground state varies only 
weakly over the whole range of pion masses and is about one half of the physical value at 
the critical point. The excited states are influenced by the thresholds and vary much more 
strongly. 

In the remainder of this subsection, we calculate the expectation values of the 2N and 
3A" kinetic energies and some properties of the 2N and 3A^ wave functions. While these 
quantities are technically not observables, they shed some light on the structure of the 
three-body states. 

In Fig. 121 we show the expectation values of the kinetic energy for the triton ground and 
first excited states and for the two-nucleon states as a function of the pion mass M^^. All 
expectation values are evaluated in the rest frame of the corresponding states. The triton 
ground state kinetic energy stays fairly constant as is varied. The kinetic energy of the 
first excited state, however, approaches the kinetic energy of the two-nucleon bound state 
(the deuteron for M.^ < M^^** and the spin-singlet deuteron for > M^"^) near the value 



8 



10.0 



> 

CD 



1.0 



0.1' \ I I I X I 

190 195 200 205 210 
[MeV] 
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deuteron (Mt^ > M^"^) indicated by the sohd hnes as function of pion mass M^^. 



°iio 

Q. 

I 

(U 

> 

o 1 

CO 

I 

C\J 

0.1 

170 180 190 200 210 

[MeV] 

FIG. 4: 2N-3N wave function overlap for the triton ground state and the first two excited states 
as function of pion mass M,r- The vertical dotted line indicates the critical pion mass MJ''**. 
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of the pion mass where the triton excited state enters from the 2N-N continuum. As a 
consequence, the third particle is essentially at rest in this region. This behavior indicates 
that the first excited state has a 2N-N cluster structure close to the 2N-N threshold. A 
similar observation applies to the second excited state which is not shown in the figure. 

The 2N-3N wave function overlap for the triton ground state and the first two excited 
states is shown in Fig. E] as function of pion mass M^. The overlap rapidly approaches zero 
as the pion mass reaches the critical value. This is a consequence of increasing the size of 
the deuteron (M^r < M^"*) and spin-singlet-deuteron (Mjr > M^"^) as the critical point is 
approached from below and above, respectively. Directly at the critical point, they are both 
infinitely large. While new shallow three-body states appear at values of that form a 
geometric series and differ by factors of Aq = 515 as M^"* is approached, the size of a given 
three-body state remains finite (cf. Fig. |2)). Consequently, the overlap between 2N and 3A^ 
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wave functions vanishes at the critical point. 

We have also calculated the probability for the D-wave and mixed symmetry states in 
the triton ground state and the first two excited states as function of pion mass M^. These 
results will be discussed below in the context of the contact EFT. 



III. CONTACT EFFECTIVE FIELD THEORY 

In principle, one could calculate all three-nucleon observables in the vicinity of the critical 
pion mass in the chiral EFT. However, the computational effort is significant, even in the 
3N system. Moreover, the calculations become increasingly difficult for the shallow excited 
states near the limit cycle. The physics near the limit cycle, however, can also be described 
in the contact EFT for large scattering length. This theory is formulated in an expansion 
around the limit of two-body bound states at threshold corresponding to vanishing inverse 
scattering lengths. It is much simpler than the chiral EFT and does not contain pion degrees 
of freedom. Three-body calculations can typically be carried out on a personal computer. 
Therefore, it is an ideal tool to calculate physical observables in the critical region where the 
scattering lengths are large. We note that the contact EFT is universal and does not have 
the explicit pion mass dependence. When the pion mass dependence of the spin-triplet and 
spin-singlet scattering lengths (cf. Fig. H} as well as one three-body observable (cf. Fig. ^ 
are taken from the chiral EFT calculation, the pion mass dependence of other three-body 
observables can be calculated. Therefore, it complements the chiral EFT study from the 
previous section. 

For practical purposes, it is convenient to write down this theory in the Lagrangian 
formalism using so-called "dibaryon" fields. In our case, we need two dibaryon fields: (i) a 
field ti with spin (isospin) 1 (0) representing two nucleons interacting in the ^Si channel (the 
deuteron) and (ii) a field Sa with spin (isospin) (1) representing two nucleons interacting 
in the ^5*0 channel [Tt^ : 

-| {slN^a.r^r.N + h.c) - [^7,^(t.a.)^(t,a,) 

+ ^ {itm)\SaTa) + h.c.) + g%S aTa)\s tU)] N , (14) 

where i,j are spin and a, b are isospin indices, M is the nucleon mass, and gt, Qs, and H are 
the bare coupling constants. The ai (r^) are Pauli matrices acting in spin (isospin) space. In 
the two-body sector, the exact solution of the field theory can be obtained analytically jsj . 
Renormalization can be implemented by adjusting the two-body coupling constants 5't(A) 
and (7s (A) as a function of the ultraviolet momentum cutoff A such that the spin-triplet 
and spin-singlet scattering lengths have the desired values. Other two-body observables are 
then independent of A and have the appropriate values up to corrections of order ro/|a| and 
To >y M\E\ where tq is the range of the interaction and E the typical energy. 

In the three-body sector, the nonperturbative solution of the field theory can be obtained 
by solving generalized Skorniakov-Ter-Martirosian integral equations including a three-body 
force numerically. These integral equations have unique solutions only in the presence of an 
ultraviolet cutoff A. The resulting predictions for three-body observables, although finite, 
depend on the cutoff and are periodic functions of ln(A) with period tt/sq where sq = 
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FIG. 5: The three-body parameter A,,, determined from matching the first excited state of the triton 
as a function of the pion mass M-,^ (circles). The dash-dotted hne is a fourth order polynomial fit 
while the vertical dotted line indicates the critical pion mass M^"^. 



1.0062378... is a transcendental number. In Ref. [13] it was shown that the quantum field 
theory could be fully renormalized to remove the residual dependence on A in the three- 
body sector by adding a three-body interaction term to the Lagrangian density in (jl4|) . The 
dependence of three-body observables on the cutoff decreases like 1/A^ if the three-body 
coupling constant has the form [s^ 



H{A) 



cos[so ln(A/A*) + arctan(so)] 
cos[so ln(A/A^,) — arctan(so)] 



(15) 



for some value of A*. With this renormalization, three-body observables have well-defined 
limits as A — ^> oo, but they depend on the parameter A,,, introduced by dimensional trans- 
mutation. Since H{A) is a periodic function of ln(A), the renormalization of the field theory 
involves an ultraviolet limit cycle. This EFT has succesfully been applied to various nuclear 
three-body observables |9]. Higher order corrections can be calculated as well 3^ 3^, 3^, 3^ 



but are suppressed by ro/\a\ and ro^/M\E\ near the exact limit cycle. (For a recent formal 
study of the corrections in repulsive partial waves, see Ref. |37||.) 

As discussed above, we take the pion mass dependence of the spin-singlet and spin-triplet 
scattering lengths from the chiral EFT calculation (cf. Fig. P). The pion mass dependence 
of the three-body parameter A* can be determined from matching the energy of the triton 
ground state or one of its excited states (cf. Fig. Since the higher-order effects are 
generally smaller for the shallower states, we match to the first excited state. The pion mass 
dependence of all other three-body observables can then be predicted to leading order in 
the power counting. In Fig. El the circles give the value of A,, obtained from the matching 
procedure as a function of the pion mass. The pion mass dependence of A,, is well described 
by a fourth order polynomial fit indicated by the dash-dotted line. For pion masses between 
190 and 210 MeV, the parameter A,, varies smoothly by about 15%. The 15% variation in 
A* corresponds to a variation of the triton ground state energy by about 25% over the same 
range of pion masses. In the previous work of Ref. the value of A^, was approximated 
by a constant since no three-body calculations with chiral potentials were available. 
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n 


^ [MeV] (chir.) 


5^"^ [MeV] (cont.) 


i?("-iVM"^ (chir.) 


B^r''/4^' (cont.) 


-1 




1.8437 X 10^ 




515.0 





3.7736 


3.5798 




515.0 


1 


6.9504 X 10^3 


6.9504 X 10^3 


542.9 


515.0 


2 


1.3287 X 10-5 


1.3495 X 10~^ 


523.1 


515.0 


3 




2.6202 X 10"^ 




515.0 


4 




5.0874 X 10-11 




515.0 


5 




9.8779 X lO-i'i 




515.0 


6 




1.9179 X 10-1^ 




515.0 



TABLE I: Binding energies of the triton excited states and their ratios at the critical pion 
mass. The second and 4th columns show the results from chiral EFT while the third and 5th 
columns show the results from contact EFT. The dashes indicate entries that have not been cal- 
culated, while the blank entries are not defined. 

We are now in the position to calculate the pion mass dependence of the other three- 
nucleon observables. The results for triton ground state and second excited state are com- 
pared to the chiral EFT results in Fig|21 The circles (ground state), squares (first excited 
state), and diamonds (second excited state) give the chiral EFT result, while the dashed 
lines are calculations in the contact EFT. We find good agreement between the chiral and 
contact EFT calculations which is generally better for the shallower states. The typical 
error is about 10% for the ground state and below 1% for the second excited state. The first 
excited state is reproduced exactly because of the matching procedure. 

The pionless theory can also be used to calculate the binding energies of the next few 
excited states in the critical region where it is very difficult to perform calculations in the 
chiral EFT. Our results for the binding energies of the first 6 excited states at the critical 
pion mass are compared to the results from chiral EFT in Table lU The second and 4th 
columns show the binding energies and ratios from chiral EFT while the third and 5th 
columns show the binding energies and ratios from contact EFT, respectively. The dashes 
indicate entries that have not been calculated, while the blank entries are not defined. From 
the ratios of the chiral EFT results it is evident that the first two excited states are still 
influenced by higher-order effects such as the finite range of the chiral potential. The exact 
ratio of 515.035... will be reached for the shallower states, but their calculation in the chiral 
EFT is computationally very expensive. The contact EFT, on the other hand, has an limit 
cycle in the ultraviolet by construction. If the theory is tuned to the critical point, the limit 
cycle is exact for all energies. The ratio in the 5th column of Table |l]is therefore 515.035... 
for all states. The contact theory is most accurate for the shallower states where the chiral 
EFT also approaches a limit cyle and becomes less accurate for the deeper states. Note also 
that this theory predicts infinitely many deeper states whose binding energies are beyond 
the range of validity of this EFT. As an example, we have shown the state with n = — 1 in 
Table HI 

The binding energies of the excited states directly at the critical point can also be obtained 
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n3 
o 



0.1: 



0.01: 



P(D) g.s. 
P(mix) g.s. 
P(D) 1st ex. 
P(mix) 1st ex. 



170 



180 



190 
[MeV] 



200 



210 



FIG. 6: Probability of the D-wave and mixed symmetry states for the triton ground state and the 
first excited state as function of pion mass M^. The vertical dotted line indicates the critical pion 

mass M^"K 



analytically from the formula [Ji 

_ (e2Vso)i-" ^(1) _ (515.035..)^-" B^'^ , (16) 

and the numerical results of the contact EFT are in good agreement with this analytical 
formula. 

If one is interested in small deviations from the exact limit cycle, for example to calculate 
observables in the critical region, precise numerical techniques are required. Results with 5 
digits of precision as in this work can be obtained in a straighforward way. Higher numerical 
accuracy requires more advanced techniques as described in Ref. where three-body 
binding energies up to 12 digits of precision have been obtained for the bosonic problem. 
Moreover, the leading dependence of the three-body energies on the physical ultraviolet 
cutoff, provided by the long-range one-pion exchange in our case, was also given in Ref. j3^: 
A renormalization group analysis suggests that the leading corrections to the three-body 
binding energies are of order v^ln(A)/A and ^ ln(A)/A2 relative to 5^"^ itself. The 
first correction proportional to ^/B2 vanishes at the critical point. The second correction 
is qualitatively consistent with the results for the chiral EFT in Table H] and predicts a 
decreased value for B-^ / B\ . However, the decrease is too strong and would lead to a ratio 
bI^^ / B^^ < 515, so other corrections must be important. 

It is interesting to compare the probability for the D-wave and mixed symmetry states 
in the chiral EFT with the expectations from the contact EFT. In Fig. we show the 
probability for the D-wave and mixed symmetry states in the triton ground state and the 
first excited state as function of pion mass for the chiral EFT. The probabilities for the 
second excited state have also been calculated but are not shown. They behave similar to 
the first excited state. The mixed symmetry probabilities in the triton excited states become 
very small near the critical pion mass. This is a consequence of the SUjA) symmetry at the 
critical point, which leads to a decoupling of the mixed symmetry state jl7l |. 

In the chiral EFT, the S'f/(4) symmetry at the critical point is exactly realized for the 
scattering lengths, but not for the effective ranges. This suggests that the mixed symmetry 
component is induced in the contact EFT by higher order terms proportional to r^V ME. 
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FIG. 7: The S = 1/2 and S = 3/2 neutron-deuteron scattering lengths (M^ < M^"*) and the 
5 = 1/2 neutron-singlet-deuteron scattering length (Mtt > M^"*) as function of pion mass M^^. 
The vertical dotted line indicates the critical pion mass M^**** while the vertical dashed lines show 
where the second excited state appears in the spectrum. 



The observed suppression of the mixed symmetry component at the critical point is indeed 
consistent with this result of the power counting. The probabilities decrease roughly by 
factors of 30, which is in the order of magnitude expected from the binding energy ratios of 
the states. The strong enhancement of the mixed symmetry state probabilities once the pion 
mass is tuned away from the critical point is generated by SU(4) breaking of the scattering 
lengths, which is a leading-order effect in the contact EFT. We observe that the mixed 
symmetry state is very large, for the excited states. Treating SU(4) breaking interactions 
as a perturbation, this behavior might be explained by the proximity of other (virtual) 
excited states, which leads to an enhancement of the perturbative contributions to the wave 
function. 

For the D-state probabilities, similar observations hold. However, since the D-state com- 
ponent is generated by subleading interactions for all pion masses, we do not observe the 
strong enhancement, once the pion mass is tuned away from the critical mass. All D-state 
probabilities decrease slowly with increasing pion mass. This indicates that the tensor force 
becomes less important at larger pion masses. This can be expected, since in the limit of very 
large pion masses, the pion becomes a heavy degree of freedom, so that the approximation 
by contact interactions becomes more and more accurate. 

The calculation of three-body scattering observables in the contact EFT is also straight- 
forward i9j]- In Fig- HI we show results for the 5 = 1/2 and 5 = 3/2 neutron-deuteron 
scattering lengths for M^^ < M^"^. Above the critical pion mass the deuteron ceases to 
exist, but the spin-singlet-deuteron becomes bound. For > M^^**, we also show the 
5 = 1/2 neutron-singlet-deuteron scattering length. The vertical dashed lines show where 
the second excited state appears in the spectrum. The critical point is indicated by the 
vertical dotted line. 

All scattering lengths diverge at the critical point. The 5 = 3/2 scattering length diverges 
because it is simply a constant times the spin-triplet scattering length: 1.179 ass^ The 
5 = 1/2 scattering lengths are very sensitive to the appearance of new excited states at 
threshold as the inverse two-body scattering lengths approach zero. They jump from — cxo 
to +00 whenever a new state appears. 
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The possible existence of a limit cycle in the pion mass region around 200 MeV must also 
be taken into account in future nuclear lattice calculations. If the pion mass is in the critical 
region, three-body scattering observables will show a strong pion mass dependence. Since, 
however, the physics of the limit cycle can be captured using EFT methods as demonstrated 
in this work, it should not render such calculations impossible. One can simply rely on the 
EFT to extrapolate through the critical region to physical pion masses. A similar strategy 
was followed in a recent fully dynamical lattice calculation of the nucleon-nucleon scattering 
lengths |3S]. 

IV. SUMMARY & CONCLUSION 

In this paper, we have performed the first study of the conjectured infrared limit cycle 
in QCD in the framework of chiral EFT. We have exemplified three parameter sets in the 
chiral EFT at NLO that lead to an infrared limit cycle for pion masses around 200 MeV 
and investigated the structure of three-body observables for one of these sets in detail. 

Using both chiral and contact EFT's, we have calculated the energies and structure 
of the triton ground state and the first two excited states around the critical pion mass. 
Furthermore, we have calculated the next four excited states and the neutron-deuteron and 
neutron-singlet-deuteron scattering lengths in the critical region. All three-body scattering 
lengths diverge at the critical point. The S = 1/2 scattering lengths are very sensitive to 
the appearance of new excited states in the critical region. 

Moreover, we have elucidated the consequences for future three-body lattice QCD cal- 
culations. On the one hand, scattering observables are very sensitive to new three-body 
states appearing at threshold in the critical region. On the other hand, the energies of the 
three-body states themselves and the triton ground state energy in particular change only 
very slowly as one passes through the critical region. This suggests that chiral EFT could 
be a reliable tool for extrapolations of lattice results from unphysically large pion masses to 
the physical value. 

The comparison of the chiral and contact EFT results also showed that an accurate 
matching of both theories is possible around the critical point. The small deviations of the 
contact EFT predictions from the full chiral result were in line with the expectations from 
the power counting of the contact EFT. It is reassuring to confirm these expectations by an 
explicit calculation. 

In summary, the main findings of Ref. [l^j have been confirmed and extended. Future 
lattice studies promise interesting insights into the rich structure of QCD. In particular, it 
would be very interesting to see signatures of the limit cycle in lattice simulations of the 
three-nucleon system with pion masses around 200 MeV j4l|. 
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